Spatial correlations of Rydberg excitations in optically driven atomic ensembles 
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We study the emergence of many-body correlations in the stationary state of continuously-driven, 
strongly-interacting dissipative system. Specifically, we examine resonant optical excitations of Ryd- 
berg states of atoms interacting via long-range dipole-dipole and van der Waals potentials employing 
exact numerical solutions of the density matrix equations and Monte-Carlo simulations. Collection 
of atoms within a blockade distance form a "superatom" that can accommodate at most one Ryd- 
berg excitation. The superatom excitation probability saturates to | for coherently driven atoms, 
but is significantly higher for incoherent driving, approaching unity as the number of atoms in- 
creases. In the steady state of uniformly-driven, extended one-dimensional system, the saturation of 
superatoms leads to quasi-crystallization of Rydberg excitations whose correlations exhibit damped 
spatial oscillations. The behavior of the system under the van der Waals interaction potential can 
be approximated by an analytically soluble model based on a "hard-rod" interatomic potential. 
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I. INTRODUCTION 

Strong, long-range dipole-dipole (DD) or van der 
Waals (vdW) interactions between atoms in highly ex- 
cited Rydberg states [l| can suppress multiple Rydberg 
excitations within a certain interaction (blockade) vol- 
ume, while enhancing the rate of single collective excita- 
tion [l-lll. The resulting dipole blockade constitutes the 
basis for a number of promising quantum information 
schemes and interesting many-body effects involving 
long-range correlations and crystallization of Rydberg ex- 
citations [13413 ■ 

Much of the research on strongly interacting Rydberg 
atoms is focused on the unitary dynamics or eigenstates 
of the many-body system. In coherently driven ensembles 
of atoms, depending on the strength and detuning of driv- 
ing lasers, different ground-state phases with crystalline 
order emerge |10l - [l5| . Yet, adiabatically attaining the 
ground state of a large system requires exceeding prepa- 
ration times [1^ during which the decoherence and 
dissipation associated with the optical excitation cannot 
be neglected. This necessitates the consideration of open 
systems and their stationary states 27[ , which is the main 
purpose of this paper. 

We study theoretically resonant optical excitations of 
Rydberg states of atoms interacting with each other via 
DD and vdW potentials. Atoms within a blockade dis- 
tance form a "superatom" which can accommodate at 
most one collective Rydberg excitation 0, H, [28l - l30| . We 
show that the steady-state excitation probability of the 
superatom saturates to ^ for coherently driven atoms, 
but can be significantly higher in the presence of strong 
dephasing which suppresses inter- atomic coherences and 
disentangles the atoms ^3d\. This case is amenable to 
an iterative Monte-Carlo sampling algorithm which can 
accurately and efficiently simulate the stationary state 
of the many-body system. We apply this algorithm to 



extended one-dimensional (ID) atomic ensembles, and 
explore both the low-density (lattice) regime [Ij, [13] 
with one or few atoms per blockade distance, and the 
high-density (continuous) regime 1233 with many atoms 
per blockade distance. For the vdW interacting atoms 
at high densities, tight packing of superatoms leads to 
quasi-crystallization of Rydberg excitations, similar to 
the predictions of an analytic model involving a "hard- 
rod" (HR) interatomic potential [13]. In contrast, the 
DD potential appears to be too "soft" for assigning a 
well-defined interaction range at high atomic densities. 

The paper is organized as follows. In Sec. |lT]we intro- 
duce the Hamiltonian and dissipative terms of the master 
equation for the density matrix of the many-body system 
and examine the properties of superatoms under coher- 
ent and incoherent driving. In Sec. lIIII we study extended 
ID systems, employing exact solution of the density ma- 
trix equations for several atoms, and Monte Carlo sim- 
ulations for hundreds of atoms. A model based on the 
HR interatomic potential is analytically solved in Sec. lIVI 
and compared with the results of numerical simulations 
for the vdW interacting atoms. Experimental considera- 
tions and conclusions are summarized in Sec. |V] 



II. THE MANY-BODY SYSTEM 

We consider an ensemble of TV atoms irradiated by 
a uniform driving field that couples near-resonantly the 
atomic ground state \g) to the highly excited Rydberg 
state \r) with Rabi frequency il, Fig. [IJa). A pair of 
atoms i and j at positions x.^ and Xj excited to states \r) 



interact either via the DD (p = 3) [3l| or vdW {p = 6) 
fs^] potential 7iA(xi — Xj) = hCp\xi — Xj\~P . In the frame 
rotating with the driving field frequency w, the system 
Hamiltonian H = Vat + Vaa is composed of the atom- 
field and atom-atom interactions. Vat = i^^rr + 
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FIG. 1. (a) Level scheme of atoms interacting with the driv- 
ing field Q on transition \g) \r) with detuning 5, while 
Vaa denotes the DD or vdW interaction between the atoms 
in Rydberg state |r) having the (population) decay rate Tr- 
(b) Atoms within the blockade distance db form a superatom 
with at most one Rydberg excitation, (c) Mean number of 
Rydberg excitations na of superatom containing A^sa < 10 
atoms within L — 0.7db (vdW interaction), obtained from the 
exact steady-state solutions of Eq. ([1} for 5 = 0, Fr- = O.lfi, 
and F2 = (filled blue diamonds) and F^ = f2 (filled red cir- 
cles). Also shown the Rydberg excitation probability (Srr) 
of superatom (empty red circles) as per Eq. (|4]). 



(7^^)] and Vaa = ^J2i<i Kr^i^i ^ ^3)Hn where 



rg 



ajj^^ = \fJ-)jj{v\ is the transition {/j, ^ v) or projection 
(/i — v) operator for atom j at position Xj, and 5 — 
uj — ujrg is the driving field detuning. The relaxation pro- 
cesses affecting the atoms include the spontaneous (ra- 
diative) decay of the excited state \r) with rate F^, and 
the (non-radiative) dephasing of atomic coherence drg 
with rate F^; the decay rate of the Rydberg state is typ- 
ically small compared to $7, while the physical origins of 
dephasing include non-radiative collisions, Doppler shifts 
or the excitation laser linewidth. The decay and dephas- 
ing Liouvillians, acting independently on each atom, are 



given, respectively, by C^p 



]. The den- 
sity matrix p of the A''-atom system obeys the master 
equation 



and Cip = F4(ct^^ - al )p(al^ - a' ) - /5] [3 



dtp^ ~-[n,p]+ Cp, 



(1) 



with Cp^j:f{Cip + Cip) 



A. Blockade distance 

For a single two-level atom, the steady-state popula- 
tion of the excited state \r) is a Lorentzian function of 



detuning S, 



{(Trr) 



(2) 



with the width W — Jrg^/'^^'^ /^rjrg + 1, whcrC Jrg = 

' 2F2 is the total (transversal) relaxation rate of the 



iF 



arg coherence [33]. For strong (fi^ > Trjrg) resonant 
{S ^ w) driving, the population saturates to ((Trr) ^ ^■ 
But given an atom in the Rydberg state |r), it will in- 
duce a level shift A — equivalent to detuning 6 — of an- 
other atom, blocking its Rydberg excitation when A > w. 
We may therefore define the blockade distance dh via 
A(db) = w, which yields 



(3) 




with p = 3 for the DD interaction and p 
interaction. 



: 6 for the vdW 



B. The superatom 

Consider iVga atoms within distance L < d^,, such 
that A(xj — Xj) 3> w for any pair of atoms i and j. 
Fig. [Ifb). We thus expect that the "superatom" can 
accommodate at most one Rydberg excitation |2^, j2^ . 
This is confirmed by our exact numerical simulations for 
-^sa < 10; starting with all the atoms in the ground 
state \g), we propagate Eq. ([T]) for time t long 
enough until the steady-state is reached. In Fig. [TJc) 
we show the resulting mean number of Rydberg excita- 
tions riR 
riR 

double (a. 



(X]^^" ^rr) within the superatom. Clearly, 
< IVTVsa, while we verify that the probabilities of 



' vrh 



triple 



etc. excitations are 
always small. 

In the absence of dephasing, Y z — ^ and 7rg = ^Fr <C 
ri, we have n^^, — \ independent of A'sa- The ground 
state \G) = 52, • ■ • , ffJVsa) of the superatom is cou- 
pled only to the single collective Rydberg excitation state 



/WZ^i |<7i,<72,-.-,'j,. 
the states ji?^"-') with higher number n > 1 of Rydberg 
excitations are shifted out of resonance by the strong 
interatomic interaction Vaa and therefore are not popu- 
lated. As the collective Rabi frequency V-^sa^^ saturates 
the transition \G) -H> the superatom ground and 

excited states acquire populations (Sgg) — {^rr} — \i 



J, while all 



where Egg = Hj^i ^Ig and J^rr 



n 



i¥^3 99 

are projectors onto the ground and single Rydberg ex- 
citation states of A^sa atoms. We note that nR slightly 
larger than ^ seen in Fig. [IJc) is due to imperfect block- 
ade of Rydberg excitation [A(L) ^ lOui] of atoms at the 
boundaries of the region of finite size \L = Q.ldh]- 

Remarkably, strong dephasing F^ increases the mean 
number of Rydberg excitations tir > ^ within the su- 
peratom. The transversal relaxation 7rg ^ destroys 
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FIG. 2. Rydberg state populations {a^j.) oi N = 3,5,7,11 
atoms (from top to bottom: red, green, blue, black) in a ID 
lattice of length L = 3db (vdW interaction) obtained from the 
exact solutions of Eq. ^ for 5 — 0, Tr — O.IO, and = 
(upper panel) and = Q (lower panel). Also shown the cor- 
responding Rydberg excitation probabilities a^^ (empty cir- 
cles) obtained from MC simulations. 



the inter- and intra-atomic coherences, causing individ- 
ual atoms to behave independently. In the basis of collec- 
tive, singly-excited states, this corresponds to a coupling 
of the symmetric state l-R^^-*) to all the non-symmetric 
states [301 ■ The superatom still contains at most one 
Rydberg excitation, T^gg + ^rr = Ij which, upon com- 
bining with a^gg + al^ = 1 , yields 



RR) 



(4) 



where {arr) is given by Eq. ([2]). For A'sa = 1 we have 
(Tirr) = (cTrr) as it should, while A'sa ^1/ (o'rr) leads to 
(Tirr) — >■ 1. InFig.[TJc) we plot {T,rr), which reproduces 
well the exact numerical solution, riR ~ (Y^rr), under the 
same conditions. 



III. EXTENDED ID SYSTEM 

We next consider the steady-state distribution of Ry- 
dberg excitations in a ID system of A' = PatL atoms of 
linear density pat- 



Simultaneously, we observe an onset of spatial oscilla- 
tions of Rydberg excitations ((J^^). Note that at higher 
atomic densities {N > 7), the oscillations are somewhat 
smoother and more pronounced in the presence of strong 
dephasing jrg ^ ^■ 



B. Monte-Carlo algorithm 

The dephasing suppresses interatomic coherences and 
disentangles the atoms, admitting only classical iV-body 
correlations. Each atom then behaves as a driven two- 
level system of Eq. ([2|) but with the detuning 6 deter- 
mined by operator Sj = X^i^j <5'rr^(^i ~ ^j) which de- 
scribes the total interaction- induced shift of level \r) for 
an atom at position Xj involving the contributions of all 
the Rydberg atoms a*^ at positions x^. 

We can now introduce an efficient procedure to sim- 
ulate the stationary distribution of Rydberg excitation 
probabilities at any atomic density pat- Our algorithm re- 
lies on iterative Monte-Carlo (MC) sampling of {a^^} 
an ensemble of N atoms, in the spirit of the Hartree-Fock 
method. We start with, e.g., all the atoms in the ground 
state, {(^ig) = IVj e [1,^^], although the resulting 
steady-state does not depend on the initial configuration. 
At every step, for each atom j, we draw a uniform ran- 
dom number s € [0, 1] and compare it with the Rydberg 
state population (al^); if s < (ct^^), we set tr^^ 1, oth- 
erwise a-j.^ — >■ 0. In turn, the thus constructed binary con- 
figuration of Rydberg excitations {al^} — > {0, 1, 0, . . .} 



determines the level shift Sj (equivalent to detuning S) 
of atom j when evaluating (a^j.)- We continuously iter- 
ate this procedure, sifting repeatedly through every atom 
in the potential generated by all the other atoms, in the 
self consistent way. The probability distribution a-j.^ of 
Rydberg excitations results from averaging over many 
10^) configurations {a^r}- We note a related algo- 
rithm employed in [s^] • In Fig. [2] we plot for various 
N and verify that, for strong dephasing {Tz = fi), the 
algorithm accurately reproduces the exact steady-state 
populations (ir^r)- the absence of dephasing (F^ = 0), 
however, the steady-state resulting from the MC simu- 
lations quantitatively differs from the exact state of the 
system. 



A. Small system 

We have performed exact numerical simulations of the 
density matrix Eqs. ^ ioi N < 11 atoms with and 
without dephasing F^. In Fig. [2] we show the station- 
ary populations (a^r) of Rydberg states of vdW interact- 
ing atoms in a lattice of length L — Sdb [open bound- 
ary conditions (OBC)]. When the interatomic distance 
exceeds {N — 3), interactions play no role and the 
population (ir^^) for each atom is given by Eq. With 
increasing atomic density, interactions progressively sup- 
press the Rydberg state populations of individual atoms. 



C. Large system 

We employ the above MC procedure to study realisti- 
cally large ID systems of up to A' ^ lO'^ atoms interacting 
via the DD and vdW interactions and subject to strong 
dephasing. In Fig. we show the spatial distribution 
of Rydberg excitation probabilities in a finite system of 
length L = 15db (OBC) for different atomic densities. In- 
teractions between the atoms suppress the Rydberg exci- 
tation probabilities (t^^, the more the higher is the atomic 
density (smaller the interatomic distance) and stronger 
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FIG. 3. (a) Probabilities al,. of Rydberg excitations in a ID 
atomic ensemble of length L — 15dh (OBC) obtained from 
MC simulations with 5 = 0, Fr = O.lfi, and Fz — Q. Trian- 
gles (connected by dotted lines) correspond to the DD inter- 
action, and circles (connected by solid lines) to the vdW in- 
teraction between the atoms. The atomic densities are, from 
top to bottom: patdb = l|, 2|, 5, 10, 20 (red, green, blue, 
brown, black), (b) The corresponding ID spatial correlations 
g^^\x) of Rydberg excitations of DD (left) and vdW (right) 
interacting atoms. In the inset of each graph, g^^\x) is plot- 
ted versus x in units of the corresponding collective blockade 
distance deb of Eq. ((5]). 



the interaction. The Rydberg excitations are then re- 
pelled to the boundaries of the system. An atom at the 
edge, having higher probability to be excited, suppresses 
the excitation of the neighboring atoms within the block- 
ade distance, beyond which another atom acquires higher 
excitation probability. Hence a Rydberg excitation den- 
sity wave develops. This behavior is significantly more 
pronounced for the vdW interaction between the atoms, 
which is stronger within the blockade distance but falls- 
ofF fast outside of it, as compared to the DD interaction, 
which is weaker at short distances and has longer tails. 
The vdW potential is thus more reminiscent of a "hard- 
rod" interaction, while the DD potential is much "softer" . 

From many configurations {a-f.^}, we can extract the 
spatial correlations of Rydberg excitations G*-^^ {x = |xi — 

Xjl) = al^air- More explicitly, we place a Rydberg atom 
((T*^ = 1) at the origin = 0, and then use the above MC 
procedure to calculate tr^^ Vx^ > 0, from which we obtain 
the normalized correlation function as g^-^^x) = al^/urr^ 
where = SjLi ^Ir the spatial average with N ^ 
1 (L >• db)- In Fig. [3Jb) we show the corresponding 
^(^^(x) for the DD and vdW interactions. In the case 
of vdW interaction, the Rydberg atom at a; = almost 



completely blocks the excitation of all the atoms within 
the blockade distance x dh- In contrast, for the DD 
potential the excitation blockage is only partial as 5^^^ [x) 
grows nearly linearly in the region of a; < db. In addition, 
for higher atomic densities Pat^b > 5, we observe damped 
spatial oscillations of g'^'^\x) with increasing amplitude 
and slowly decreasing period A close to db- 



D. Collective blockade distance 

In a dense system, p^tdh ^ !> an atom in the Rydberg 
state |r) blocks the excitation of other atoms within a 
certain distance which, due to collective effects, is some- 
what smaller than the blockade distance db for a pair 
of atoms, Eq. ([3]). To estimate the collective block- 
ade distance deb, note that the collective Rabi frequency 
f^cb ~ ViVcbf^, and thereby the excitation linewidth 
Wch — Sr^cb VTrg/Tr f^r A'cb = Patdcb atoms is enhanced 
by a factor of \/Nch- Substituting iVcb into the definition 



of deb = ^Cp/wch yields 



deb — 



2/(2p+l) 



(Patdb)l/(2f+l) ■ 



(5) 



with p — 2> for DD and p — Q for vdW interactions. In 
turn, the number of atoms within the collective blockade 
distance is iVeb = {pM^"''^^"^^^ ■ 

In the insets of Fig. |31[b) we plot the correlation 
functions g'^'^\x) with x rescaled by the corresponding 
(density-dependent) collective blockade distance deb- For 
vdW interaction, all the curves for different atomic densi- 
ties then exhibit the same oscillation period A ~ 1.75dcb- 

The average (background) density of Rydberg excita- 
tions Pvdw = PatO'rr Can also be deduced from the col- 
lective excitation picture. For the average probability of 
the Rydberg state of vdW interacting atoms we obtain 



(iVsa 
\12/13 



l)(a. 



1' 



(6) 



where iVsa ~ 1.83{p^tdh) ' ~ 2iVeb is the effective 
number of atoms per superatom [see Fig. Ufa)]. Note 
that superatoms overlap. 



IV. HARD-ROD POTENTIAL MODEL 

We have seen above that the vdW interaction between 
the atoms in Rydberg states results in almost complete 
blockade of simultaneous excitation of two or more atoms 
within a distance close to db. This suggests an analogy 
between the vdW interaction potential and a HR poten- 
tial of a range dr ~ db [21]. The model of incoherently 
driven atoms interacting via the HR potential admits an 
analytic solution presented below. Our aim is quanti- 
tative comparison of its predictions with the result of 
numerical simulations for vdW interacting atoms at high 
densities. 
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FIG. 4. (a) Average density pvdw = patO'rr of Rydberg ex- 
citations of vdW interacting atoms vs. the atom density pat 
obtained from MC simulations (circles) and Eq. (|6]) (solid 
red line). Also shown the average excitation density pun of 
Eq. (|10|l for the HR potential of range dr = db (dotted black 
line), (b) Oscillation period A and (c) decay length ^ of spa- 
tial correlations of excitations g^'^\x) vs. the atom density 
Pat- For vdW interaction potential, the data points (red cir- 
cles with uncertainty) are extracted from the MC simulations; 
for HR potential, A and ^ are obtained from the solution of 
Eq. (|15|l . The effective range di^'^^ (blue stars; right vertical 
axes) of the HR potential is obtained by equating A (b) and 
^ (c) for the vdW and HR interactions. 



A. Rate equations model 

In we have introduced a rate equations treatment 
of the ID system with complete blockade of Rydberg 
excitations of neighboring atoms. We now extend this 
formalism to the longer-range interatomic interactions. 
To this end, we consider N atoms on a lattice, with 
one atom per site, and assume that the atoms excited 
to the Rydberg state \r) interact via the HR poten- 
tial of range di-: A(xi — Xj) = CooB(|xi — Xj| — dj) 
with Coo S> w and Q{x) being the step function. An 
atom in state |r) blocks the excitation of = patdr 
neighboring atoms on both sides along the chain. For 
each atom, we then have two incoherent processes: a 
pump from \g) to |r) with rate P, conditioned upon 
the absence of Rydberg excitations within the range dj.. 
Lip = riixi-ajj Kd, (^rr " 1); and a deexcitation 



from |r) to \g) with rate D, h 
the two rates k = -£ = , 

Li i — \cr^, 

with (5 = as 



d — 



Dal,. 



The ratio of 



is obtained from Eqs. ([2]) 



101' 



2 TrTrg 



(7) 



The density operator of the system obeys the equation 
of motion 



5*p = E 



{L^L^,p}+2Ll,pL 



(8) 

After sufficient relaxation time, the density matrix at- 
tains a classical form p = > 
where p{{pj}) is the probability of the TV-atom configu- 
ration {p,j} e {g,r)^ . In the steady state of the system, 
we have the detailed balance relations 



M-M' 



(9) 



where M = '^rr llMil) is the total number of 

excitations in {pj}. States with the same number of ex- 
citations have equal weight and the partition function is 
given by Zn = Y.m=oHM,N)k^' , where A(M,iV) = 

(^~^M^~^'') is number of possible arrangements of 
M excitations on a lattice of N sites, with any two excita- 
tions separated by at least iV^ sites. The mean number of 
excitation is given by M = X]a/=o ^(-^' N)k^ while 
the average (background) density is phr = M/N. In a 
large system N ^ N,-, A(A/, 7V)k*^ is a highly peaked 
function of M; finding its maximum at M^ax — M we 
obtain the density 



Phr = 



1 WjP) 
d, 1 + W{p) ' 



(10) 



where j3 — kN^. and W{l3) is the Lambert function de- 
fined via We^ — /?. For /3 ^ 1, the density approaches 
Phr — >■ We note that Eq. (ITOl) is also obtained 

in the classical model of hard rods [33 where dr is the 
rod length and /3 is a free parameter. Our derivation, 
however, yields 



/3 = Kpatdr, 



(11) 



which is uniquely determined through the system param- 
eters. 

Next to the boundary at a; = 0, the density of excita- 
tions reads 



Phr(x) ^ , ^e(a;/d..~fc) 



dr 



k=0 



!i^{x/dr - kf 
fc! 



(12) 

and the spatial correlations are given by g'^'^^x) — 
Pbr{x — dr)/pYm- At distances x > d^, the correlations 
have the form of decaying spatial oscillations 



g(2)(a;) ~ 1 + Acos(27ra;/A + 0)e-^/«, 
where the oscillation period and decay length, 

X^dr-, 

C = drhi"^ 



WiP) sin(&) 



(13) 

(14a) 
(14b) 
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are determined by the solution of the transcendental 
equation for b, 

ln(^(/3))+M.(/^)-ln(-i^)+^.0. (15) 

B. Comparison of the vdW and HR potentials 

To quantify the suitability of the HR model [21| for the 
description of the system under the vdW potential, we 
now compare the mean densities of excitations Pvdw.HR, 
as well as the oscillation periods A and decay lengths ^ 
of spatial correlations g^^''{x). 

Let us first equate the range of the HD potential to 
the blockade distance of the vdW potential, dr = dh- 
The average densities of excitations are then remarkably 
close, Pvdw — Phr, especially at high atomic densities 
Patdh ^ 1, see Fig. Hl^a). With increasing the atomic 
density, however, the oscillation period A of the density 
wave decreases faster for the vdW potential. Fig. HJ^b) ; 
apparently the superatoms having soft boundaries can 
pack closer to each other. The same softness of the vdW 
potential, compared to the HR potential, can explain 
shorter correlation lengths ^, see Fig.[3]^c). 

Next, we may assume equal oscillation periods A for 
both potentials and then deduce the corresponding range 
di^^ of the HR potential, which turns out to be close to db 
but slowly decreasing with increasing the atomic density 
Pat, Fig. mjb). Alternatively, we equate the correlation 
lengths ^ and find somewhat smaller corresponding range 
of the HR potential, dP ~ 0.8db, again slowly decreasing 
with increasing the atomic density pat, Fig.UJc). 

Finally, we have performed MC simulations for ^ 1 
strongly-driven atoms interacting via the HR potential. 
In Fig. [5] we compare the the correlation functions g'^^^x) 
obtained from our simulations and the corresponding an- 
alytic solutions with the results for the vdW interacting 
atoms under the otherwise identical conditions. We ob- 
serve qualitatively similar behavior when the range d^ of 
the HR potential is equal to the blockade distance d^ of 
the vdW potential. Alternatively [2l|, we may associate 
the range of the HR potential with the collective block- 
ade distance deb in which case /3 = KA^cb = ^(Patdb)^^^^^. 
As seen in the insets of Figs. [SJ the correlation functions 
for vdW and HR potentials are now considerably differ- 
ent, the corresponding oscillation periods being an ap- 
proximately constant A ~ l.TSdcb for the vdW potential, 
and a slowly varying A ~ 1.3 — 1.25dcb for the HR poten- 
tial. 



V. CONCLUDING REMARKS 

The system studied in this paper corresponds to an 
ensemble of cold alkali atoms excited to the strongly- 
interacting Rydberg states |r) with principal quantum 
number n 50 — 100. The resonant atomic excitation 
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FIG. 5. Correlations g^'^\x) of atomic excitations for vdW 
(circles) and HR (triangles) potentials obtained from MC sim- 
ulations with atom densities /Oatdb ~ 30 (top) and patdb — 50 
(bottom). In the insets, g^^^{x) are plotted versus x in units 
of the collective blockade distance deb- Atomic parameters 
are as in Fig. [3] 



is effected by either direct one-photon (UV) transition 
Ig) — >■ |r) or two-photon transition via non-resonant in- 
termediate state |25| . Typical values for the driving field 
Rabi frequency fi ^ 10^ Hz, together with the relaxation 
rates F^ ^ O.lfi and jrg — 2r2, lead to the blockade dis- 
tance in the range of db ~ 5 — 10 pm. The low-density 
ensemble of equidistant atoms is trapped in a ID optical 
lattice. At higher densities, regular arrangement of atoms 
in a lattice would play a minor role and our results should 
hold also in the continuous gas of atoms confined in an 
elongated trap with transverse dimension much smaller 
than the blockade distance. 

To summarize, we have shown that the steady-state 
probabilities of Rydberg excitations in resonantly-driven 
atomic ensembles can exhibit damped spatial oscillations 
and quasi-crystallization due to the strong interactions 
between the atomic Rydberg states. The inverse wave- 
length and correlation length of these oscillations grow 
with the probability of collective Rydberg excitations of 
superatoms which increases with the atomic density and 
thereby the cooperativity of the excitations. An analytic 
model based on a "hard-rod" interatomic potential can 
serve as a guide to understanding the properties of the 
system under the van der Waals potential. 

After sudden switching-off of the driving field, the Ry- 
dberg quasi-crystal can survive for tens or hundreds of 
microseconds, it can be detected in situ by spatially- 
resolved Rydberg state ionization [2^ or high-resolution 
fluorescence imaging [2^ . 
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